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Abstract 

We consider the asymptotic behaviour of positive solutions u of the conformal scalar 
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curvature equation, Au + ' ~ ' u n ~ 2 = 0, in the neighbourhood of isolated singularities 
in the standard Euclidean ball. Although asymptotic radial symmetry for such solutions 
was proved some time ago, we present a much simpler and more geometric derivation 
of this fact. We also discuss a refinement, showing that any such solution is asymptotic to 
one of the deformed radial singular solutions. Finally we give some applications of these 
refined asymptotics, first to computing the global Pohozaev invariants of solutions on the 
sphere with isolated singularities, and then to the regularity of the moduli space of all such 
solutions. 

1 Introduction 

The problem we consider in this paper is to derive asymptotics for positive solutions of the 
conformally invariant scmilinear elliptic equation 

. n(n — 2) 

Au+— =0 (1) 

4 

which are defined in the punctured ball B n \ {0}, and which are singular at the origin. It is 
well-known that a solution u of this equation corresponds to a conformally flat metric 

g = U n - 2 

which has constant scalar curvature R g = n(n — 1). Equation (|l|) is a special case of the more 
general equation relating the scalar curvatures of any two conformally related metrics. More 

4 

specifically, if g and g' — u n ~ 2 g are any two such metrics, with corresponding scalar curvature 
functions R g and R g i, respectively, then 

n — 2 n — 2 n+2 , 

4(n — 1) 4(n — 1) 
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L 9 = A 9 - TU—- H R 9 ( 3 ) 



The linear part of the operator on the left 

2 

4(n -^T) ' 

is known as the conformal Laplacian associated with the metric g. It satisfies the conformal 
equivariance property that 

Lg(u<j>) = U^Lg,(f> (4) 

for any cj) G C°° and any two metrics g and g' related as above. 

Asymptotics of solutions of equation ([!]) were proved by L.A. Caffarelli, B. Gidas and J. 
Spruck § using a rather complicated version of Alexandrov reflection. Their result states that 
an arbitrary positive solution u of ([!]) satisfies 

u{x) = \xfPv e (r log \x\ + T)(l + o(l)), (5) 

where u e (x) — \x\^~ v £ (— log |x|) is one of a one-parameter family of radial solutions of ([[]). We 
call the v e and the u s Fowler solutions and discuss them at greater length in the next section 
(Fowler solutions are also called Delaunay type solutions in many works which we will refer to). 
In recent work C.C. Chen and C.S. Lin have generalized this result and slightly simplified 
their argument, although their proof is still quite technical. This asymptotics theorem is related 
to a simpler result due to P. Aviles [|l| concerning the subcritical equation, where the exponent 
2i| is replaced by a smaller number p in the range [—35, ;rr§)- P- Aviles shows that an arbitrary 
positive solution of the subcritical equation with an isolated singularity at the origin converges, 
in the sense of (J5J) , to a radial solution of the same equation. 

The present work is a combination of two separate projects which are sufficiently closely 
related that it has seemed natural to publish them together. The first, undertaken several years 
ago by the first and fourth authors (K & S), is to supply an alternate argument, more geometric 
and much simpler than that of 0, to establish (||). This also gives a slightly stronger estimate, 
improving the o(l) remainder term to 0(|a;| Q ) for some a > 0. The second project, quite recently 
undertaken by the second and third authors (M & P) and inspired by their recent construction 
0, concerns a refinement of these asymptotics, in some sense obtaining a second term in the 
expansion, along with various applications for these refined asymptotics. More precisely, it is 
possible to improve (j^) using, instead of v e , a family of deformations of the Fowler solutions, 
parametrized by vectors a S K", which we denote v e ^ a and which are defined by 

u e ,„(t,0) = \9 - ae-'l^v^t + log\9 - ae-'l), 
where 9 = x/\x\. For all o£R™ and TeRwe define 

2 — n 

U e ,a,T(x) = \x\~V S: a(-\og\x\ +T,x/\x\). 

The main result is 

Theorem 1 Let u be an arbitrary weak solution of the equation (jl|) which is positive on the ball 
B n (0) C R n and in C co (B n (Q) \ {0}). Then either u € C°°(B n (0)), or else there exists some 
choice of parameters e € (0, £0] f £ o = (( n ~ 2)/n) ( ™~ 2 - ) / 4 / ), a € l n and some constants T € R and 
a > 1 such that 

u(x) = \x\^(v e , a (-log\x\ +T,x/\x\) + 0(\x\ a )), as \x\ -> 0. 

Alternatively, 

u(x) = u e . a .T{x) + 0(\xf ) as |ie|-»0 
for some constant [3 > ^-j^ . In particular, there are upper and lower bounds for u of the form 

2 — 71 2— Tl 71 + 2 

Cijxj - ~ <u< C2\x\^~ , withO < Ci < C2, and the corresponding metric g — u n ~ 2 5 is complete 
at 0. 
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Solutions of equation (||) on the sphere S n , singular along a closed set A, giving conformal 
factors relating the standard metric to a new metric of scalar curvature n(n — 1) complete on 
S n \ A, were first considered in the work |^0) of the fourth author and S.T. Yau. For solutions to 
exist it is necessary that the Hausdorff dimension of A be less than or equal to (n — 2)/2. Under 
mild geometric hypotheses, these solutions are shown there to extend to global weak solutions on 
S n . Around the same time, solutions with isolated singularities (and with some more complicated 
singular sets arising essentially as limit sets of Kleinian groups) on S n were constructed by one of 
us [ p~7| . Since then, solutions singular along an arbitrary disjoint collection of submanifolds, and 
giving metrics of positive scalar curvature, have been constructed in a succession of works, |l3| ], 
PI , § , |§ ■ (We do not discuss the by now extensive set of results on the analogous problem for 
metrics of negative scalar curvature.) 

Emerging from || and j^] (as well as |l^]) is the lesson that solutions with isolated singularities 
are far more rigid objects than solutions with higher dimensional singular sets. This is apparent 
if one fixes the singular set A and studies the moduli space Ma of all solutions singular along 
A. We shall henceforth restrict attention to the case where A is a finite disjoint collection of 
submanifolds, of (possibly varying) dimensions less than or equal to (n — 2)/2. When at least one 
component of A is of positive dimension, Ma is infinite dimensional and is modelled locally on a 
Banach manifold. Closely related to Theorem [l] is a result on refined asymptotics, analogous to 
the one here but for solutions singular along higher dimensional submanifolds, given in [jij. One 
sees from these asymptotics that solution metrics need not have bounded curvature, and that the 
correct way to parametrize the moduli space Ma then is via an asymptotic Dirichlet problem. 
On the other hand, when A is a set of k points on the sphere, then Ma is finite dimensional. It 
is proved in [[llj that in this case, Ma is even a locally real analytic set with virtual dimension 
k. In H it is shown that for generic configurations A, Ma attains this dimension. The sharper 
asymptotics proved here have applications to the study of Ma when A is finite. In particular, 
corresponding to any conformal Killing vector field on the sphere one can define a Pohozaev 
invariant of the (global) solution, depending on the parameters e and a associated to each singular 
point. Using the refined asymptotics one can evaluate these invariants explicitly; one obtains a 
finite number of real analytic equations satisfied by the set of parameters associated to each 
solution. While these equations are not enough to determine the moduli space completely, they 
do lead to new information about it. This is done in §6. There are two other applications of the 
refined asymptotics considered here. Then, in §7, we discuss the nondegeneracy of the moduli 
space Ma near its ends. 

For any k > 0, Theorem [j] gives the existence of a 'parameter map' $ on the space of solutions 
defined on all of K™, with exactly k singularities; <!> associates to a solution its set of singular 
points pi , . . . , pk and the associated parameters Tj , aj and Sj at each pj . $ is defined more 
carefully in (|2l]) below. This map should be quite useful in the global study of the moduli spaces 
Ma when A is finite. It seems quite likely that <E> is injective, but the proof of this fact may be 
subtle. We leave this as an interesting open problem. 

We conclude this introduction by remarking on the relationship of the results here with what 
is known for noncompact constant (nonzero) mean curvature surfaces ScK 3 which are (Alexan- 
drov) embedded and have k ends. Although seemingly quite different, this extrinsic geometric 
problem bears many similarities to the intrinsic conformal scalar curvature equation. Indeed, all 
the results known in one setting seem to have direct analogues in the other setting, although the 
proofs are sometimes different. Any end of such a surface £ is asymptotic to one of the Delaunay 
surfaces (these are a family of rotationally symmetric CMC surfaces in R 3 ); this follows from 
work of the first author, Kusner and Solomon, Q, and Meeks, [Q. The deformations of the 
radial solutions used here, and the corresponding refined asymptotics, corresponds in the CMC 
setting to translations and rotations of the asymptotic axes. More analytically, the difference 
is that the Jacobi fields corresponding to these deformations, in the CMC context, are at most 
linearly growing, while in the scalar curvature setting they are exponentially growing or decaying. 
Amongst the various applications of the refined asymptotics we develop here, the global balanc- 
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ing formula of §6 corresponds to the computationally more simple balancing formula in the CMC 
context discussed, for example, in [Q. The nondegeneracy argument of §7 has a direct analogue 
in the CMC setting, and to our knowledge, is new there as well. 



2 Preliminaries 

In this section we define and discuss in some detail the (n + l)-parameter family of functions i> Ej£1 
which arise in the statement of Theorem |l| We also study the linearization of ([j]) about any one 
of the the Fowler solutions v e , since this will be a crucial tool in the later analysis. Parts of this 
material can also be found in 1 7 and [[flj . 



In the ensuing discussion it will be convenient to use the conformal equivariance of the problem 
to express (|l|) relative to different background metrics in the flat conformal class. The three most 
natural choices are the standard Euclidean metric, the cylindrical metric dt 2 + dd 2 on K x S*™ -1 , 
and the standard spherical metric on S n . Equation ([I]) is already expressed relative to the 
Euclidean metric. Relative to the cylindrical metric this equation becomes 



d 2 v . (n-2) 2 n(n-2) »+ 



, ,Agv-- — -t^ v + ^~a — -v*=*=0. (6) 



The transformation between ([[]) and (0) can also be seen directly by setting 

2 — n 

u(x) = \x\^~v(— log |x|,x/|a;|) 

and choosing a new independent variable t = — \og\x\. We shall not need the explicit form of 
this equation relative to the spherical metric. 

At each stage below we shall use whichever of these ambient metrics is most convenient for 
the computations at hand. 

2.1 Fowler solutions 

We now define the special family of solutions mentioned in the introduction and develop their 
relevant properties. We start by considering the radial solutions of (|l|) on R" which are singular 
at the origin. For this problem it is simplest to use the cylindrical metric and equation <^). 
The first result characterizes the global solutions of (^). 

Proposition 1 Let v be any positive solution of defined on the whole cylinder M x S n ~ 1 . 
Then v is independent of the spherical variable 9, hence depends only on the variable t. Moreover 
we always have < v(t) < 1 for all t € R. 

The proof of this uses an Alexandrov reflection argument, and is contained in [|| and also in |jq| . 

Hence if v is a global solution, then it satisfies the autonomous ODE obtained by dropping 
the 9 differentiation in (||). We analyze this by converting it to a system of first order equations 
on the phase plane. Namely, setting w = dtv, the equation becomes 

d d (n-2) 2 n(n-2) 2+2 

— V = w, —w = V v n ~ 2 . 

dt dt 4 4 

This system is Hamiltonian, with corresponding energy function 

tt, \ 2 {n-2) 2 2 (n-2) 2 _2n_ . . 

H(v, w) =w 2 - —v 2 + —v "-2 . (7) 

v ' 4 4 

In particular, if v(t) is a solution of then the path parametrized by (v(t), v'(t)) is contained 
within a level set of H, hence H(v,v') does not depend on t. There is a homoclinic trajectory 
lying in the level set {H = 0} fl {v > 0}, and a one-parameter family of closed level sets contained 
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in the bounded set O = {H < 0} n {v > 0}. All other trajectories and level sets pass into the 
region where v < 0, so we do not consider them here. The compact level sets in O correspond 
to a family of periodic solutions of (H) which we shall denote v £ (t). The parameter e denotes 
the minimum value attained by the solution, and we will often refer to it as the necksize of the 

4 

solution. It varies in the interval (0, £o], where £q ~ 2 = IL -^ L - The level set of H corresponding to 
v e is 

H^H{e) = ^^{e^-s 2 ). (8) 

Before continuing with the rest of our treatment of this special family of solutions we digress 
slightly to discuss the energy of arbitrary (nonnegative) solutions of (JsJ) . Thus let v be any such 
solution. Multiply this equation by dtv; a small calculation yields 

U(d t v) 2 +divg(d t vVov) - klVH 2 - { -^ld t (v 2 ) + = 0. 

2 zoo 

We can rephrase this by noting that it is equivalent to the vanishing of the divergence of the 
vector field 

W = (±{d t vf ~ \\V e v\ 2 - ^^v 2 + ^l v ^ dt + (d t v)V e vV 9 . 

Hence the integral of (W, dt) on the sphere t — T (with respect to the volume form for the metric 
g) gives a number which, by the divergence theorem, does not depend on T. Thus we obtain an 
invariant of the solution v, which we shall call the 'radial Pohozaev invariant' of v and denote 
Vra^v). It is a special case of the more general Pohozaev invariants we shall discuss later. Note 
that the radial Pohozaev invariant of the Fowler solution v £ is simply ^uj n _iH(e), where w„_i is 
the volume of S"™" 1 and H(e) is the Hamiltonian energy ((7J). 

We now return to the discussion of the special family of solutions. There are actually two 
parameters for these periodic solutions. The first is this necksize parameter e, while the second 
corresponds to the value of the solution at t — and we will denote by T e the period of v e . 
We normalize the functions v e by assuming that v E (0) — min?; e . Hence the complete family 
of positive radial solutions is v £ (t + T) for T € K and e <G (0, £o]. The two extremes of this 
family are when v is the constant solution v = eo (geometrically, this corresponds to the correct 
magnification of the cylinder which has scalar curvature n(n— 1)), and the limits v E (t) — ► and 
v e (t + \T e ) — > (coshi)^ - as e — + 0. This last explicit solution corresponds to the conformal 
factor (coshi) -2 transforming the cylinder to the punctured sphere S n \ {p, —p}. 

It is known that T e , the period of v e , is monotone in e, converging to 2tt / \fn — 2 as s — > eo an d 

increasing to oo as e — * 0. Geometrically, however, the metrics g e = Vs~ 2 (dt 2 + d0 2 ) converge to 
a bead of spheres arranged along a fixed axis. This family of radial solutions interpolates between 
the cylinder and this singular limit. 

We may also transform these back to solutions of (0), to obtain the solutions 

e^ T u £ (e- T x) = \xfrv e {-\og\x\ +T) 

on M™ \ {0}. 

As indicated earlier, we shall consider, in addition, an n parameter family of deformations of 
these radial solutions. Geometrically these may be understood as follows. Using the Euclidean 
background metric, and starting with a radial solution u, we first take the Kelvin transform 
u(x) = \x\ 2 ~ n u(x/\x\ 2 ). It is not difficult to check that u is also a solution of Then translate 
u by some vector a £ R™, and finally take the Kelvin transform once again to obtain u e , a . Note 
that this is the pullback by the composition of three conformal transformations, first inversion 
in the unit sphere {|x| = 1}, then Euclidean translation, which is a parabolic transformation 
fixing infinity, and then inversion once again. Each of these elementary transformations takes the 
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space of solutions of ([j]) to itself. Using instead the spherical metric as the background, if u is 
transformed to a solution on S n which is singular at two antipodal points {p, —p}, then these steps 
correspond to first reflecting across the equator determined by p and —p, then applying a parabolic 
transformation F which fixes the point p and carries another point q to p, and finally reflecting 
across the equator again. Finally, note that relative to the cylindrical background metric, the final 
and initial transformations correspond even more simply to the reflection t — > — t. The parabolic 
translation is more complicated to describe on the cylinder. 

Now let us perform these steps explicitly on u £ (x) = \x\^~ v £ (— log |x|). After the first inver- 
sion, we obtain |x| ~~ ~ v £ (\og \x\). Translation carries this to \x — a\~^~ v £ (Yog \x — a\). Performing 
the final inversion yields at last 

X 2 — n X 2— n 2 — n 

u e , a (x) = | — - a\ — v £ (\og\ — - a\) = \x\ — \6 ~ a\x\\~ v e (- log |x| + log|0 - a|x||), 
p| \x\ 

where 9 = x/\x\. We can also define the associated functions 

« e ,«(t,0) = Iff-ae-^WeCt + logie-oe-*!) 

on the cylinder; these are solutions of (||). Notice that these functions are regular except at 
t = log \a\,8 = a/\a\, and in particular are smooth for t sufficiently large. 

We conclude this section by calculating the expansion of u £j£I near |x| = 0. To do this we first 
observe that 

l^-aMl 2 ^ = 1 + ^—^a ■ x + 0(\x\ 2 ). 
'|x| 2 

Similarly we calculate 

log It— r - a|x|| = —a ■ x+ Of|x| 2 ), 
|x| 

and so 

v e (- log |x| - a ■ x + 0(\x\ 2 )) = v e (- log |x|) - v' e {- log |x|)(a ■ x) + 0(|x| 2 ). 
Putting these altogether we finally obtain 

2-r> / n — 2 \ 

u e , a (x) = \x\— f u e (-log|x|) + a-x(-4(-log|x|) + — — ^(--ogM)) +0{\x\ 2 )\ 

2 — n Tl - — 2 6 — n 

= u e (x) + \x\ — (a ■ x)(-v' e + —j-**) + 0(\x\ — ). (9) 

In particular, \v e>a (t,6) — v e (t)\ < Ce~ t . 

We have now obtained the families of solutions u £M (x) and v £M (t, 9) for the problems ([!]) and 
(^), respectively. For all Tel, changing v E (t) into v e {t + T) in the construction above leads 
to families of solutions u ei a,T(x) and v e ^T(t,0) for the problems (0) and (^), respectively. We 
note again that the parameters a £ M. n and T £ K correspond to explicit geometric (conformal) 
motions; only the parameter e does not arise from an extrinsic motion. For simplicity, when 
a = 0, we shall denote these functions simply by u e and v s . 

2.2 The linearized equation 

We shall now consider the linearization 

d 2 (n-2) 2 n(n + 2) ^ 
^ = ^2+ A « ^ + 4 v - ( 10 ) 

of the operator in (||) around one of the radial Fowler solutions v £ . Our primary interest is in the 
mapping properties of C £ , which we shall review from JTTJ and |l. 
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This operator has periodic coefficients, hence may be studied by classical Flocquet - or Bloch 
wave - theoretic methods, as in 111]], but also by separation of variables and elementary ODE 
methods, as in |)). We shall need to refer to results derived from both of these methods. 

Let {Xj, Xj(8)} be the eigendata of A S n-i. We use the convention that the eigenvalues are 
listed with multiplicity, so that Ao = 0, Ai = . . . = A„ = n— 1, A n +i = 2rt, etc. Then C £ decouples 
into infinitely many ordinary differential operators 

d 2 ,n(n + 2) « (n - 2) 2 , 
C ^ = dT 2 + V " -^f— "^)- 

When j > n + 1, the term of order zero in each of the C £ j is negative because Xj > 2n and 
v £ < 1; hence for these values of j, £ EJ satisfies the maximum principle. The same conclusion 
holds for C e j, j = 1, . . . , n, because conjugating by an appropriate power of v £ yields an operator 
also with negative term of order zero. It follows from these facts that the L 2 spectrum of —C e j 
is contained in (0, oo) when j > 0. On the other hand, is contained in the essential spectrum 
of — £ e .o- In particular, C does not have closed range on L 2 (R x S" -1 ; dtd9), but the difficulty is 
localized to the ground eigenspace of the Laplacian on the cross-section. 

The failure of C £ to have closed range is caused by its Jacobi fields, i.e. by the solutions of 
C £ tp = 0. To study these, it suffices to consider the solutions of the induced problems C e jt/jj — 0. 
It can be proved, cf. pi, that for each j > there are two normalized linearly independent 
solutions ipfj and an associated constant 7 £J - such that < e TlE - jt for all t. The solutions 

ifj^ (t), ip~ (t) for C e fl are bounded and linearly growing, respectively, and so it is natural to 
define 7 £ ,o = 0. We immediately deduce the following 

Proposition 2 Suppose that £ e ?p — on K + x S 71 ^ 1 and ip — 0(e~ 7t ) as t — > oo, for some 
7 > 0. Then ip(t, 9) — Y^jLja ^Ej{^)Xj(^)> where jo is the first integer such that 7 Jo > 7. 

Somewhat remarkably, we may deduce the explicit forms of these solutions when j = 0, . . . , n. 
This rests on the fact that whenever we have a one-parameter family of solutions of the nonlinear 
equation (j^), then the differential of this family gives a solution of the linearized equation. In the 
previous subsection we found several different one-parameter families of solutions. The simplest 
are the families 

T — >v £ (t + T), and e — ► v e (t). 

Differentiation with respect to either of these two parameters yields the solutions corresponding 
to the cross-sectional eigenvalue A = 0. 

<o(*) = ^v e (t + T)= v' e (t), ip- (t) = ^v s (t) = v s (t)- (11) 

We may also differentiate v e ^ a with respect to ay, using the expansion (^), to get 

d n — 9 

^(t)xi(9) = ^v e , a (t, 9) =e-\- v' £ (t) + —v e (t))xj(0); (12) 

these give the exponentially decreasing solutions of of C £ jip = 0, j = l,...,n. Since these 
functions are all equal to one another, it is often simpler to denote any one of them by just V^i- 
In any case, using them we may restate (^|) as 

u e , a (x) = \x\ ^ (v e (- log |x|) + N (E "SXiWtfA- l0 S I^D) + °^ x \ 2 ) j 

= \x\^ (v e (- log \x\) + (a ■ x)^+x(- log \x\) + 0{\x\ 2 )) . (13) 

Notice that ipti = ■ ■ ■ = tytn ■ Although they are not required later, we also note that the 
exponentially increasing solutions r/>~ -(t) of £ £ ,j, j = 1, . . . , n, are obtained from tptji^) simply 
by replacing t by —t. 
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It follows from ( |12[ ) that the 'indicial root' 7 e .j = 1, j = 1, . . . , n. It is somewhat unex- 
pected that for j < n, j s j is independent of e; one can show however that all other 7 ej depend 

nontrivially on e, and in fact tend to (- — — — — h Aj) 1 ^ 2 as e — > ||. 

The numbers {±7 ej } are analogous to the indicial roots of a Fuchsian operator, and they 
influence the mapping properties in exactly the same way. To state these mapping properties 
most suitably for our purposes we define weighted Holder spaces on the (half-) cylinder 

Vt G R, C^ a {[t , +oo) x S"- 1 ) ={w= e^w : w G C k > a ([t , +oo) x 

If T is any space of functions defined on the half-cylinder [to , +oo) x S n ~ 1 , we let Tx> denote the 
subspace of functions of T vanishing at t — to for some fixed to- We cannot choose to arbitrarily 
because it is necessary for the following result that tf>j(to) ^ for j < n; this condition is always 
fulfilled for j = 1, . . . , n but fails for j = when t € {{l/2)T £ :feZ}),T £ being the period of 
v £ . We also define 

E efi = Span{V+ (t)} C C 2 "([io, +oo) xr 1 ) 

and 

E £ . x = Span{^+.(t) Xj -(0) : J = 1, . . . ,n} C C\°([t , +oo) x S"- 1 ). 
Using this notation, the result we need is 

Proposition 3 Assume that to ^ (£/2)T £ , I G Z. Recalling that 7 e .o = 0, 7^1 = . . . = 7 ei „ = 1 

and 7e, n +i > 1, ifO < 7 < 1, iften 

£ e : [C% Q ([t ,^) x S»-i)e£ £i0 ] c ^C°.«([t 0)OO ) x 
is a surjective Fredholm mapping. In particular, there exists a bounded right inverse 
G £ ,o : C°'"([io,oc) x S"- 1 ) — > [C!^([*o,oo) x 5"" 1 ) 8 ^,0]^ 

SO i/iai CsGe.O = I- If 1 < 7 < 7e,n+l> 

£ £ : [C% Q ([i ,oo) x r^B.offi^j -^C°^([t ,oo) x S"- 1 ) 
zs afoo a surjective Fredholm mapping, and there exists a bounded right inverse 

g eA : C% Q ([i ,oo) x — ► [C% Q ([i ,oo) x £ e , © E e ,i] c 

with C £ Q £ .i = I . 

One may find solutions of C £ w = f for / G C_'"([to) +00) x by considering the induced 

differential operators £ £ j on each eigenspace of Ag separately. When j > n, solutions are easily 
obtained using the fact that C £ .j satisfies the maximum principle. This is done explicitly in Q. 
To obtain solutions when j < n it is easiest to use the Fourier-Laplace methods of jll| . Although 
these methods are L 2 -based, it is not hard to check that w G C_'"([io, +00) x S"^ 1 ) if / G 
Cj~([to, +00) x S*™ -1 ). The necessity of adding E £i o and E e> i to the domains to obtain surjectivity 
for the homogeneous Dirichlet problem is a simple form of the linear regularity theorem of fil|| , 
but also follows from easy ODE arguments. 

It is important in this result that the weight factor —7 does not equal one of the ±7 e ,j- In 
fact, on any one of the spaces C±™ B . ([to, +00) x S*™ -1 ) the operator C £ does not have closed range. 

This proposition may be used not only to find solutions of the equation C £ u = /, but also to 
obtain decay properties of solutions which are already given. 

Corollary 1 Suppose that fi and 7 are weight parameters and C £ u — f where u G C 2 .'^ t (R + x 
S"- 1 ) and f G C"-"(M+ x S"- 1 ). J/0 < /i < 7 < 1, then u G C 2 '"(R+ x S*"- 1 ). IfO < fj, < 1< 
7 < 7 £ ,„+i ; then u G C^(R+ xS"" 1 )®^. 
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The proof of this corollary is straightforward. For example, if0</x<7<l, fix some to > 
satisfying the assumptions of Proposition ||, then v + ap£ = Ge,of <= [C^'™([to> +00) x S*™ -1 ) 
-£'e,o]'D is also a solution of £ e u = /. Hence h = u — v is a homogeneous solution, £ e /i = 0, which 
is also exponentially decaying as t — > +00. This means that ft g Clf ([t , +00) x 5™ 1 ), and so 
we conclude that u has the stated decay. The other cases are proved similarly. 



3 Upper and lower bounds for solutions 

In this section we derive upper and lower bounds on solutions defined in the punctured unit ball. 
We further show that any such solution has a nonzero radial Pohozaev invariant "P ra . Specifically 
we prove the following. 

Theorem 2 Assume that u is a nonnegative smooth solution of ^) defined in the punctured unit 
ball B n \{0}. Either u extends as a smooth solution to the ball, or there exist positive constants 
C\, C2 such that 

cm^' 2 <«(*) <c 2 \x\ (2 - n ^ 2 . 

Furthermore, the radial Pohozaev invariant of u is nonzero. 

As stated, the constants C\ and C2 depend on the solution u. We will see in this section that the 
upper bound C2 follows from a more precise universal upper bound, while the discussion of the 
previous section, and in particular the fact that the infimum of v £ tends to zero with e, shows 
that the lower bound depends in a more delicate way on the solution. It seems to be unknown 
how one should express the lower bound in a universal form. A universal upper bound for global 
weak solutions was given in Q, and was later recorded and used in [jl6|. A local version for 
smooth solutions was obtained by the fourth author, and presented in his course at the Courant 
Institute in 1989-90. This result was written up and applied by Y.-Y. Li ||. 

The main subtlety in obtaining an upper bound is the existence of spherical solutions on K™. 

Definition 1 For any A > and Xq € R n , the spherical solution of dilation A and center xq is 
given by 



u\,x (x) = \x - xqI" 2 " (cosh(-log|x -x \+ log A)))' 



2A 



Any one of these functions gives a conformal factor transforming the flat metric on K n to the 
pullback of the standard metric on the sphere by a conformal transformation. The danger is that 
a general solution may be well approximated by one of these, in particular by a strongly dilated 
one, in some neighbourhood. This corresponds to the phenomenon of bubbling. 

We present here a local version of the upper bound for solutions defined in a region which 

extend to be weak supersolutions across the complement. A positive function u in L" o ~ 2 [B n ) is 
a weak supersolution of (|l]) if it satisfies 

n(n — 2) n+2 

Au + — -u~ < 

4 

in the distributional sense on B n . We will be especially interested in supersolutions which are 
smooth solutions on an open subset Q. of B n . The complement A = B n \il is then relatively 
closed in B n . Note that such a function u is superharmonic, and thus can be redefined on a set 
of measure zero so as to be upper semicontinuous. As such, the restriction of u to any compact 
subset of B n achieves its infimum. 

We remark that in the main application, u will be a weak solution. The reason for considering 
supersolutions is that it was shown in ( |2fj|| , Theorem 5.1) that any solution u which defines a 
complete metric on a region VL in the sphere extends as a weak supersolution to the full sphere. 
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n + 2 

Moreover, any solution u on O which lies in L™ ~ 2 (B n ) and tends locally uniformly to infinity on 
approach to A extends as a weak supersolution on B n . This can be seen by observing that for 
any large constant L, the function ul = min{u, L} is a weak solution of the inequality 

n(n — 2) n+2 
Awl < ^ u^X{u<l} 

where \e denotes the characteristic function of a set E. A simple application of the dominated 
convergence theorem then shows that u defines a weak supersolution. 

The proof of the upper bound estimate for a solution u in B n which is singular only at the 
origin will require that u extends to a weak solution of ([!]) on all of B n . This follows from a more 
general result. In fact, any solution of ([j]) defined outside a 'sufficiently thin' set automatically 
extends as a weak solution. To make this precise, we let A C B" be a relatively closed set of 
Lebesgue measure zero. We will call A thin if there is a sequence fa of smooth functions on B n 
with values in [0, 1] such that fa = in a neighborhood of A, lim fa = 1 on B n \ A, and 

lim / {\Afa\ {n+2 ^ A + \Vfa\^l 2 )dx = Q 

for any r < 1. One sees that a smooth submanifold of dimension less than is thin by choosing 
fa to be a function of distance to the submanifold; in particular, a point is thin in any dimension. 



Lemma 1 Suppose u is a solution of ^) defined on B n \A where A is a relatively closed thin set. 
Then u lies locally in L^ n+2 ^ ("~ 2 ) (B n ) , and defines a weak solution of ^) on B n . 

Tl+2 

To prove this we first show that u is locally in L™^(B n ). Let ( be a smooth nonnegative 
function with compact support in B n \ A, and multiply (|l|) by £(™+ 2 )/ 2 and integrate by parts to 
obtain 

c („ + 2)/2 u(n+ 2)/(„-2) dx = _ j uA(C («+2)/2 )da ._ 

This implies the bound 

j (C ( "- 2)/2 w) (n+2)/( "- 2) dx < c j C ( "- 2)/2 w(C|AC| + |VC| 2 )dx. 

An application of Holder's inequality then implies 

J (C ( "~ 2)/2 w) (ll+2)/( "- 2) dx < c J (C|ACI + |vc| 2 ) (n+2)/4 . 

Now we choose £ = tpfa where ip is a function which is equal to one on B(0, r) for some r < 1, 
and equal to zero outside B(0,ri) for some r± G (r, 1). We then have by an easy estimate 

(C|AC| + |vC| 2 ) (n+2)/4 < c(l + |A^|(™+ 2 '/ 4 + |V^| (Tl+2)/2 ) 

with a constant c depending on r, r\ . Letting i tend to infinity we then conclude that u is in 
(B(0,r)) for any r < 1. 

To complete the proof of the lemma, we now let £ be any chosen smooth compactly supported 
function in B n , and we multiply ([!]) by C,fa and integrate by parts to obtain 

/ («A(C&) + ^^c<M ( " +2)/(n - 2) )^ = o. 

We expand 

A(Cfa) = (A(fa) + 2(VC, Vfa) + faA(. 
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By Holder's inequality 



u(A((f>i)\ < \\u(\\n n+2 )/(n-2)\\A(cj)i)\\ L (_ n+ 2)/A( spt Q, 

and this term tends to zero as i tends to infinity. Similarly, 

U<VC,V&>| < ||w|VC||| L (n + 2)/(„- 2 )||V(^)|L(n +3 ,/4 (sptC) 

which also goes to zero. Therefore we may apply the bounded convergence theorem to let i tend 
to infinity and conclude 

(„A(C) + n(n ~ 2) CM ("+2)/("-2) )dx = o. 
B n 2 

This completes the proof of Lemma [l]. 

In the following proposition we let d(x,A) denote the distance from x to A for x E B n . 

Proposition 4 Let u be a positive supersolution of (jl|) which is a smooth solution on an open 
set f2 C B n . Then there exists a constant c depending only on n such that 



for allxE £2n .8(0, 1/2). 



u(x) < cd(x.A) 2 inf u 

\3B(0,3/4) 



Let A c denote the neighborhood of radius e about A. For x E .8(0,5/8), set U e (x) = 
d e (x)^~ u(x) where d e (x) = min{c?(a;, A e ), | — |x|}. In order to prove the proposition it suf- 
fices to show that 



sup U e < c inf u) (14) 

5(0,5/8) \dB(0,3/4) J 

for some constant c not depending on e. 

The first step in proving this is to show that given any constants R and 5 E (0, 1/2], there 
is a constant c\ depending only on n, i?, and (5 such that if M e = sup B( -Q 5 / 8 ) U e > c\, then the 
solution u "differs by at most 5" from a spherical piece "of size i?" which is concentrated near 
any chosen maximum point xq E B(0, 5/8)\A e of U e . Precisely we mean that, if Xo is a maximum 
point of Ue, and we set A = u(xo) 2 ^ 2 ~ n \ then the rescaled function w\(x) = \( n ~ 2 ^ 2 u(xo + Xx) 
satisfies 

11^,^ - W\\\c*(B(0,R)) < S (15) 

where u^. Vo is the standard spherical solution defined above with \y$\ < C2 and X/c% < [i < C2 for 
a constant C2 depending only on n. To prove this assertion, we use an indirect blow-up argument. 
Assume that there is a sequence Ui with corresponding values Mi^ e tending to infinity. Let x% be 
a maximum point of Ui_ t , and let Wi(x) — A,-™ 2 ^ 2 Ui(xi + Ajx) where A^ = Ui(xi) 2 ^ 2 ^ n \ Since 
Xi is a maximum point of Ui ye , we have Ui ye (x) < U^ e (xi) for x E B(xi,ri) where r$ = ^d^xi). 
Since di t£ (x) > \di >e (xi) for x E B(xj,r,) we have Ui{x) < 2( n ~ 2 ^ 2 Ui(xi) for x E B(xi,ri). We 
then see that Wi is a solution of (|l|) which is regular and bounded by 2("~ 2 )< /2 in B(0, Ri) where 
Ri = Ti/Xi = \M 2 ^ n 2 \ Since Ri tends to infinity, we see in a standard way that a subsequence 
of the Wi converges in C 2 norm on compact subsets to a global positive solution w of ([!]). By the 
theorem of || this must be a standard spherical solution. 

We now argue by contradiction. If the result were not true, there would exist a sequence Ui 
with corresponding values Mi ;€ tending to infinity such that (|l5| ) does not hold for any /i > and 
any yo £ lR n - This clearly contradicts the previous analysis. 
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Once ( |l5| ) is proved, the bounds on y and /i follow from the conditions w\ < 2^™ 2 ^ 2 in 
B(0,R'), where i?' = ±4(x )/A > ±c 2/( "~ 2) and w A (0) = 1. Indeed, increasing c\ if necessary, 
we may always assume that c\ > 2 n ~ 2 so that R 1 > ^c 2 ^™ -2 '' > 2. It follows from w\ < 2( n ~ 2 ^ 2 
and S < 1/2 that U Mo O) < 2("- 2 )/ 2 + 6 < 2"/ 2 for all y S 5(0,2) C 5(0, i?') and it follows from 
w\(0) — 1 that u^ : y (Q) > 1 — J > 1/2. Now assume that y G -B(0, 2), we may then take y = yo 
in the inequality u fl ^ yo (y) < 2™/ 2 to obtain 2/i < 2™/(™~ 2 ) which gives a bound on \i. Finally, if 
i/o ^ B(0, 2), there exists i/i € £?(0, 2) such that | j/i — yo\ 2 = \yo\ 2 — 1 and it then follows from 
the inequality u fl ^ yg (yi) < 2™/ 2 that 

2 - 2 /(n- 2 ) M < 1 + M 2 |yo|2 _ M2 

and from the inequality Ufj, jVO (0) > 1/2 that 

l + M 2 M 2 <2"/(- 2 V. 

These inequalities imply first that 2~"/( n - 2 ) < ^ < 2"/(™- 2 ) and then |y | < 2 n /("~ 2 ). 
Notice that the conclusion ( jl5| ) can be replaced by the simpler inequality 

||"1,0 - W X \\c 2 (B(0,R)) < S ( 16 ) 

if we allow ourselves to shift the center point xq to a new point x± which is within distance 2ci\ 
where C2 = ci(n) is the constant above, and which is a local maximum point of u. Thus we now 
are taking w\(x) = \( n ~ 2 ^ 2 u(x\ + \x) with A = u(xi) 2 ^ 2 ~ n K To see this we observe from the 
blow-up argument above that Wi will have a absolute maximum point near yo for i sufficiently 
large, since \y^\ is bounded and since \x is bounded from below and from above. Shifting center 
to this point, and rescaling so that the maximum value is one then produces ui t o as the limit 
of the sequence. Note that the distance 2c2A^ is an arbitrarily small constant times r% since 
\i = 2M^ Ti, and thus X\ can be used in place of xq above and we can assume, increasing 
Ci if necessary, that 2c2A^ < 1/16. 

To complete the proof of ([ll]), we first note that inf as(o,3/4) u is bounded above since one can 
see easily from the fact that u is a supersolution (as in the proof of Lemma [l]) that u is locally 
bounded in £( ,l + 2 )/("~ 2 ) (£?"). Thus if the left hand side is bounded above by the constant C\ 
(for a fixed chosen R and S) of the previous paragraph, then we are done. Therefore we may 
assume that this is not the case, and that u satisfies the above bounds with 5 chosen as small 
as desired, and R chosen as large as we like. In this case, we are going to prove that ( |l4| ) still 
holds. We now introduce the function v(t,6) on the cylinder M x S*™ -1 which corresponds to 
the function w — w\ with A = u(xi) 2 ^ 2 ^ n \ Thus we define v(t,9) = |a;|( n_2 ^ 2 u'(a;), where 
t = — log | a; | and 9 = x/\x\. It is important to note in the following argument that because the 
point Xi lies in the ball of radius 5/8 + 1/16 = 11/16 about the origin, its distance to dB(0, 3/4) 
is at least 1/16; upon rescaling, we see that the ball B(0, y^A -1 ) lies in the domain of w, and 
corresponds to a portion of B(0, 3/4) in the domain of u. From (Jlq) , we see that v is close in 
C 2 norm on [— logi?, oo) x 5 n_1 to the function vo(t, 9) = (cosht)^")/ 2 . We now fix R = e 2 so 
that -logi? = -2. It follows that d t v(-l,9) > for all 9 e S'™" 1 . Then apply the Alexandrov 
technique to v on the region [— log( ^A" 1 ), oo) x S n ~ x , reflecting across the spheres {£i} x S n_1 , 
starting with t\ very positive, and continuing as far as possible. Because dtv{—\,9) > 0, this 
procedure must end before t\ reaches —1. We will show below that because v is a supersolution, 
and the reflected function v* is a solution, there can be no interior contact point with v* < v. 
Since v is a regular solution near t = t% for t\ > —1, the Hopf boundary point lemma implies 
that dtv(ti,9) < 0. Thus the reflection argument must end because a contact between v* and v 
occurs on the boundary. In order for this to happen we must have 

uuX-log^A^M) : 9 e S 11 - 1 } < snp{v(t, 9) : t > log^A" 1 ) - 2, 9 G S^ 1 }. 
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Recalling the definitions of v and A, we see that 

mf{v(-\og(—\- 1 ) 7 9):6eS n - 1 } = 16 {2 - n)/2 inf u > 16 (2 ~ n)/2 inf u 

16 dB(x lt l/16) 3B(0,3/4) 

where the inequality holds because u is superharmonic. Finally, for 5 small and fixed, we have 
v[t, 9) < 2vo(t) for t > 0, and therefore we have 

inf u < c sup v < cA (,l_2)/2 = cm(ii) -1 

95(0,3/4) [log((LA-i)-2,oo) 

where we have used the definition of vq and A. This implies (|l4|) since C/ e is bounded above by a 
constant times u. 

In this argument we needed to know that if v* < v in a connected open set, then either v* < v 
or u* = v. To see this, observe that L(v - v*) < ( w *)(«+2)/(«-2) _ v (n+2)/(n-2) < q wea ki y . The 
standard mean value inequality then gives us the desired conclusion. This completes the proof of 
the proposition. 

For the proof of the lower bound, we require the following result. 

Lemma 2 If u is a positive solution of (Qj defined in B n \{0} with lim x _>o \x\^ n ^ 2 ^ 2 u{x) = 0, 
then u extends as a smooth solution to all of B n . 

To prove this, first note that by Lemma |l| u extends as a weak solution to B n . If we could show 
that u € L p for some p > Inj (n — 2), then standard results would imply first that u is bounded 
near the origin, and then by linear elliptic theory that it extends smoothly across 0. Let v(t, 9) be 
the corresponding cylindrical solution, and observe that the hypothesis implies that v(t, 9) tends 
to uniformly as t tends to oo. Therefore (|l|) implies that Aw > [3v for some [3 > and t > to 
for sufficiently large to. Now consider the function w = ce _v ^' + ee^* where c is chosen large 
enough that ce~^ ta > v(to,9) for all 9 6 5 n_1 . Note that w is a solution of Aw — (3w. We will 
use the maximum principle on [£o)^1 x S n _1 for any sufficiently large T to obtain the necessary 
decay estimate on v. If T is chosen sufficiently positive, then we have v(T,9) < w(T,8) since 
v tends to zero as t tends to oo while w is growing at an exponential rate. By the maximum 
principle it follows that for all to < t < T, all 9, and any chosen e > we have 

v{t,6) < ce'^ + ee^. 

Since T is arbitarily large, this inequality holds for all t > to, and since c is independent of e, we 
may let e go to zero. We have shown that e^ t v(t,9) < c for t > to- Writing this in terms of 
u we have shown that u(x) < c\x\ q for q = (2 — n)/2 + ^/]3. This implies that u £ L p for some 
p > 2n/(n — 2) as required. This completes the proof of Lemma |2[ 

We now are in a position to prove Theorem 0. The upper bound follows from Lemma [l] and 
Proposition ^, as already noted. To prove the lower bound, it is most convenient to use the 
cylindrical background, and so consider the function v(t,9). We first observe that the upper 
bound v < C2 implies a Harnack inequality of the form 

sup{u(t, 9): T < t < T + T } < c mi{v{t, 9): T < t < T + T a } (17) 

for all T > 2, and where the constant c depends on To and C%, but not on T. Define v(t) = 
J s „_! v(t, 9) d9. From the Harnack inequality we see that v and v are bounded in ratio. Suppose 
then that the singularity at is not removable. Applying Lemma ^ we see that either the lower 
bound holds or else there is a sequence U tending to 00 such that v(ti) tends to 0. 

Using (O), the Harnack inequality, standard elliptic estimates as well as the fact that v is a 
solution of (Q), we can estimate the radial Pohozaev invariant "P ra (i>), as defined in §2 (cf. also 
§6), at any one of the U. Using the fact that it is independent of i, we see that the invariant 
must vanish. Now set Wi(t,9) — v(t + ti,9)/v(ti). Some subsequence of the Wi converges in C 2 



13 



on compact subsets of R x S 71 ' 1 to a positive solution of Lw = 0. The associated function h(x) 
defined on R n \ {0} is then a positive harmonic function, and can therefore be written a|:c| n ~ 2 + b 
for some a, b > 0. This is the same as the condition that w(t, 9) = ae( n ~ 2 ^ 2t + be^^ n ^ 2 ^ 2t . Since 
w has a critical point at t = 0, it follows that a = b, and both are positive. Then, a direct 
calculation shows that 

.lim / (l(d t w l ) 2 -l\W e w^~^^w 2 +m^ i -^^ W r)de 

l— » + 00 J gn-1 2 2 o 8 

= [ -(d t wf - ( " ~ 2)2 w 2 d9 = ~to n (n ~ 2)2 ab + 

But this is a contradiction, for this is also a limit of rcscalings of the Pohozaev invariant for v, 
and hence must be zero. Indeed 

= v{ur 2 VUv) = J sn j\^) 2 - - ( -^^w 2 + v{t i )^^^ W f")d9 ? 

for i sufficiently large. This establishes the lower bound. 

To show that V Ia (v) ^ 0, we observe that by the upper and lower bounds, we can choose a 
sequence t\ tending to oo so that the corresponding translated solutions t — ► v(t + 1^, 6) converge 
in C 2 norm on compact subsets of R x S*™ -1 to a solution satisfying the same bounds and defined 
on all of 1 x S*™ -1 . By Proposition |l| such a solution must be independent of 6, and hence must 
have nonzero radial Pohozaev invariant. Since this invariant must be the same as V Ta {v) because 
of the C 2 convergence, we have shown V Ta ,{v) ^ as desired. The proof of Theorem || is complete. 

4 Simple convergence to a radial solution 

Our aim in this section is to prove 

Proposition 5 Assume that u is a nonnegative smooth solution of ^) defined in the punctured 
unit ball B n \{0} with a nonremovable singularity at the origin. For t > and 9 G S*™ -1 , let 
us define v(t,9) so that u(x) — \x\~ ~v (— log |a;|, Then, there exists a Fowler parameter 
£ € (0, so], constants T £ [0, T e ) and C > 0, and some exponent a > such that 

\v(t,9)-v e (t + T)\<Ce- at , for t > 0. 



We first show that any sequence of translates of v has a subsequence converging to one of the v £ , 
then that any angular derivative of v converges to zero. After this, a somewhat delicate rescaling 
argument due originally to Leon Simon, in a different context, gives the final convergence. 

By the results of the last section we know that < C\ < v(t, 9) < C2 for all t > 0. By standard 
elliptic estimates, we also get the uniform boundedness of any derivative \dfd7v(t, 9)\ < Cj n for 
all t > 0. 

Let {tj} be any sequence of numbers converging to 00, and define Vj{t,9) — vit + Tj,9). 
Then Vj is defined on [— tj,oo) x S n ~ 1 and satisfies (|6|) there. Using the uniform bounds on any 
derivative of Vj, we may choose a subsequence of the Vj converging in the C°° topology on any 
compact subset of R x S 11 ^ 1 . The limit function, v^, still satisfies (g), is nonvanishing because of 
the lower bound for v and is defined on the whole cylinder. By Proposition ^, the only functions 
with these properties are the translated Fowler solutions; hence we deduce that for some e and 
T, Voo {t,6) =v e {t + T). 

The radial Pohozaev invariant of v equals that of any one of the Vj , hence also equals that of 
v s (t + T). Hence the limiting necksize e is independent of the original sequence of numbers Tj 
and of the subsequence. 
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The fact that any sequence of translates of v has a subsequence converging to a ^-independent 
solution is very strong. It implies immediately that any angular derivative dgv tends to zero 
uniformly. For if this were false, then there would be a sequence of points (Tj, 9j) with tj — ► oo 
for which \dgv(Tj,9j)\ > C > 0. But then, translating back by tj and rotating 9j to some fixed 
point #o G S"" 1 to get a new sequence of solutions Vj, we can again extract some subsequence 
converging to a radial function. But this contradicts the positive lower bound on \d$Vj(0, 9 )\. 

Next, let X be any infinitesimal rotation on S'™ -1 , which we denote for simplicity by dg; 
applying it to (^) and using that it commutes with Ag, we see that dgv is a Jacobi field for 
the Jacobi operator at v, i.e. the linearization of @ at v. By the discussion above, we know 
that dgv = o(l), but since we know little about this Jacobi operator, we cannot deduce better 
decay directly. However, for some Tj — > 00 consider the corresponding sequence of translates 
dgVj and let Aj = sup t>0 \dgVj |. If this supremum is attained at some point (sj,9j) then Sj 
must stay bounded. Otherwise we could translate back further by Sj to obtain, after passing to 
a subsequence, a Jacobi field <f> for the Jacobi operator at v s (t + T), for some T, defined on all 
of K. By construction, <f> is bounded and not identically zero, but on the other hand, since it 
clearly has no zero eigencomponent relative to Ag, the results above show that it must increase 
exponentially either as t — > +00 or as t — > —00. This contradiction shows that Sj stays bounded. 
Now fix a positive integer N such that NT e > supsj and define In = {0 < t < NT £ }: here, of 
course, T e is the period of v e . Then Aj = sup Jjv |c?0i>j*|, and we obtain again, after passing to a 
subsequence, a Jacobi field <p for the Jacobi operator at some v e (t + T) which is bounded for t > 
and attains its supremum of 1 in In- 

By the results on this Jacobi operator stated earlier, <p must decay at least like e _t . Actually, 
something slightly stronger is true. Although the limit <j> is not necessarily unique, there exists a 
constant c, independent of all other choices, such that \4>(t)\ < ce _t for t > 0. To see this, write 
(j> = E"=i Cjifr+j + <j>" = <j>' + 0", where 0" e Span{fcj(t)x*(0) : j > n + 1} = E" . Since \<j>\ < 1 
for t > 0, the Cj, j — 1, . . . , n, arc absolutely bounded because they are given by the (normalized) 
inner products of <j> with Xj(@)- Hence | Y^j=i Cj^tjl — ce ~* w ith c = c(n). Next, we claim that 
10" I < ce _t for some c independent of all other choices, and that this is a consequence of the fact 
that 4>" is absolutely bounded for t > and tfi" £ E". We argue by contradiction and assume 
that there exists a sequence of Jacobi fields <j>'( G E" such that A'[ = sup t>0 e*|^'(t, 9)\ tends to 
infinity. If this supremum is achieved at a point (tj, 6i), then define 4>i(t, 9) = e li A~ x (\>"{t + ti, 6); 
these are Jacobi fields on [— fj,oo) which are bounded by e - ', and which take the value 1 at 
(0,9i). The ti clearly cannot be bounded, because the Jacobi equation is uniformly elliptic on 
any finite piece of the cylinder. Passing to a subsequence, if necessary, we obtain a Jacobi field <fi 
defined and satisfying \<f>\ < e~* on the whole cylinder. Since we also know that <j> £ E" , it must 
grow faster than e - ' either at +00 or —00, which is a contradiction. This proves the claim. 

At this stage, it is now possible to show that dgv decays exponentially. Although we do not, 
strictly speaking, need this result, we sketch the proof anyway as a warm-up to the arguments 
later. Let J N = {NT e < t < 2NT e }, and denote by v T (t,9) the translate v(t + t,9). We claim 
that if N is chosen large enough then 

Vt > 0, sup \dgv T \ < i sup \d g v T \. 

Jn ^ In 

As usual, this is proved by contradiction. If it were to fail, there would exist some sequence 
Tj — > 00 such that Bj > \Aj 1 where 

Aj = sup \dgvj\, Bj = sup \dgvj\. 

In Jn 

However, since A~ 1 dgVj converges to the Jacobi field <j> uniformly on compact sets, we see that 
dgVj = Aj4> + o{Aj) on I N U J N . Since \4>\ < Ce"*, we then get that Bj < CAje~ NT ' + \A 3 for 
j sufficiently large. However, as we have seen, the constant C is bounded independently of other 
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choices since <ft is normalized to be less than 1 on I N , so if N is sufficiently large, Ce~ NT ' < j, 
and we obtain a contradiction. This proves the assertion. It is now trivial to deduce that dgv 
decays at some exponential rate. 

We may use the sequence of translates Vj to obtain a Jacobi field in another way. Suppose 
that we have chosen a subsequence converging to some v £ (t + T). Define 

Wj(t, 9) — va (t, 9) — v E (t + T), ctj = sup \wj\, and (f>j—aT 1 Wj. 

In 

Then it is easy to check that <fij converges to a solution <j> of the Jacobi operator at v £ (t + T). This 
Jacobi field is bounded on t > 0. To see this we examine its eigencomponents with respect to 
Ag. First consider the sum of eigencomponents over all nonzero eigenvalues, which we call <fi. To 
show that <f> is bounded, and hence exponentially decaying, it is sufficient to show that dg<j) is also 
bounded for t > 0. We may assume that <j> is not identically zero, for otherwise this case is trivial. 
This function arises as the limit of a subsequence of the sequence aj 1 dgVj. If we knew that ctj 
were commensurate with Aj, the supremum of dgVj on In, i.e. if C\Aj < aj < C^Aj for some 
constants C\,Ci > 0, then we could appeal to our earlier argument. However, this must be the 
case, for otherwise a~ x dgVj would either blow up or tend to zero on In, whereas we know that it 
converges to the nontrivial function dg(f>. As for the remaining eigencomponent, with eigenvalue 
zero, we know that 

( j ) (t,6)=a + TP+ Q + a-ib- + 4>, 

where (j> decreases exponentially. We claim that a~ = 0. Intuitively this is true because this 
coefficient corresponds to an infinitesimal change of Fowler parameter. Rigorously, simply write 
Vj = v s + a.j4> + o(aj) on In, say; if a - were not to vanish, then the radial Pohozaev invariants 
of the two sides would differ, which is impossible. Thus we have shown that (f> must be bounded 
for t > 0. 

We are finally in a position to show that the difference between v(t,9) and some v £ (t + T) 
converges to zero uniformly. The subtlety here arises because the displacement T cannot be 
detected by any of the Pohozaev integrals. We use again the same sort of argument as above, 
which is due to L. Simon. Define v T (t, 9) = v(t + t, 9) and set w T (t, 9) = v T (t, 9) — v e (t). (Since 
we do not know the correct translation parameter T beforehand, we choose one arbitrarily, say 
T = 0.) The idea is to establish an improvement of approximation in the following sense. Use 
the interval In as before and fix a constant B > 0, and let t)(t) = sup Jjv \w T \. Then we claim : 

If t is sufficiently large and t](t) sufficiently small, then there exists an s with \s\ < Bt](t) 
such that jj(t + NT e + s) < ±r/(r). 

To prove this, suppose that, for fixed values of N and B, it fails. Then there exists some 
sequence Tj — ► oo such that r\j = i](rj) — > 0, but such that for any s with \s\ < Br\j we have 
rj(Tj + NT £ + s) > \r\j- Let <f>j = v^ x w Ti . We have shown already that <f>j converges in C°° 
on compact sets to a Jacobi field (f> which is bounded for t > 0. By construction, \<f>\ > 1/2 on 
NT e <t< 2NT £ , so ^ 0. Expanding </> into eigenfunctions again, we get </> = a + ipf + <f>, where 
(j> is exponentially decreasing. Note that a + , which we relabel simply a, is uniformly bounded, 
independently of the sequence, because \(f>\ < 1 on < t < NT £ . We assume that the constant 
B has been chosen larger than this upper bound. Because ip^~ corresponds to infinitesimal 
translation, this argument indicates that we should adjust v £ by some translation, and to first 
order, this adjustment should just be Sj = —r\ja; by our choice of B, this is less than Br\j. Thus 
we estimate 

w t 3 + Sj (t, 9) = v(t + Tj - rjja, 9) - v £ (t) = w T] (t, 9) - arijip+ Q (t) + o(rjj). 
In particular, 

w Tj+Sj = ri 3 4> + o{rij). 

But the function <ft decays at a fixed exponential rate and its supremum on < t < NT £ is 
bounded by one, so that if N is chosen large enough (but again, independent of our original 
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sequence Tj, etc.) then 

sup \w Tj+Sj+NTe \ = sup \w Tj+aj \=7]j sup \4>\ +o(r)j) < -j-. 

0<t<NT s NT e <t<2NT e NT E <t<2NT s 4 

But this contradicts the assumption that 77(7^ + NT e + s) > ^i](Tj), and thus we have proved our 
claim. 

It remains only to show how this approximation improvement property yields the correct 
simple asymptotics. In fact, we must prove that u)„ — > for some fixed translation a. We 
may assume, by starting at some sufficiently large t, that 77(0) is sufficiently small and that 
Br)(0) < ^NT e . Let To = and so be determined by the first step of this argument. Define, for 
any integer j > 0, 

3-1 

Oj = S i} Tj = Tj-l + 8j-i + NT £ . 

i=0 

Then by iteration, r/(Tj) < 2~- 7 77(0), \sj\ < 2~ 3 ~ 1 NT e , and hence in particular a — lirnoj is well 
defined and less than NT e . To show that a is the correct translation parameter, for any t > 0, 
let [t] denote the reduction mod NT £ , so that t — jNT e + [t] for some j > 0. Then 

w a (t,9) = v{t + a,9) - v e (t) = (v{t + aj,9)-v e (t)) + (v(t + a, 9) - v(t + a j: 9)) 

= w Tj ([t],0)+O(2- j ), 

where we have used Taylor's theorem and the uniform boundedness of dtv and dtv £ . But using 
the bound on r}{rj), we finally conclude that 

\w*{t, 9)\ < C2~\ or equivalently \w a {t, 9)\ < C'e^^H *. 
This completes our proof of Proposition |H[ 

5 Refined asymptotics 

In this last step of the proof of Theorem |l|, we improve the asymptotics one step further. We 
show that for some a £ MP, u has an expansion of the form 

u(x) = \x\ t 2 (v e (- log \x\ + T) + (a ■ x)ip+i(- log |x| + T) + 0(\xf)) , (18) 

where (3 = min{2, 7^}. 

On the other hand, we also know that the deformed Fowler solution u e<a has an expansion of this 
exact form, (j^) ; comparing these two expansions completes the proof of the main theorem. 
Using the simple asymptotics already established, write 

1*0) = |a;|^%(-log|a;|) = \x\^ (v e (- log \x\ + T)+w(- log |x|)) 

where w(t) € C^'"(R + x S" 1 ^ 1 ) for some 7 > 0. The function v(t) satisfies (g), which we write as 
N(v) = 0. Expanding this in a Taylor series about v e (t + T) gives 

n(n-2){, ."+2 n + 2 \ _ s 

C s w = — [{v e + w)— -vr 2 ^r 2 w) = Q{w). 19 

4 \ n — I J 

It is straightforward to check that if w S C^'"(IR + x S*™ -1 ) and v £ + w > 0, then Q(w) € 
C%(R+ xS"- 1 ). 

First assume that < 7 < \. Then Corollary |l| gives that w <E C^ 7 (K + x S n ~ r ), and so 
Q(w) er^l+xS"- 1 ). Continuing on, we deduce after finitely many steps that w £ Cj", (R + x 
gn-ij for some y g ^/2, 1). Applying Corollary^ again gives w S C^(M+ x S n ^ r )®E £A , where 
(3 = min{27',7 2 }). 

The optimal f3 we could expect is min{2,72}. We have proved the expansion (|l8|). 
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6 The global balancing formula 



In this section we shall present an application of the refined asymptotics theorem. We consider 
solutions of the singular Yamabe problem with discrete singular set A in the standard conformal 
class [go] on S n . More specifically, let A = {pi, . ..,Pk} C S n be arbitrary, k > 2. The solutions 

4 

we consider are functions u > on S n \ A such that u^^go is a complete metric of constant 
positive scalar curvature n(n — 1) on S" \ A. More analytically, these functions satisfy the special 
case of (^), 

(n-2) 2 (n-2) 2 ™+2 , s 

A s ^u- V - > u+ y 4 ' u— = 0, (20) 

with m > on S n \ A and u singular at the points of A. The 'unmarked moduli space' M. k 
consists of all such solutions with k singular points, but with A allowed to vary. (One may also 
consider the 'marked moduli space' A4\, consisting of all such solutions with fixed singular set A.) 
The most basic question about these moduli spaces, whether they are nonempty, was answered 
originally by the fourth author ||l7| some time ago; an alternate construction was recently given 
in |9[|. 

The refined asymptotics theorem implies that to any u € M k one can associate a set of 
parameters: the singular points pj, and then at each pj the Fowler and translation parameters 
Sj G (0, £o] and dj £ K n . This yields a map 

<S>:M k ^X k = ^[jS^UyA,, x (0,e ] fe x (M") fc j /£*, (21) 

where is the diagonal of the product of the i th factor and the j th factor and E^ is the 
permutation group on k letters, which acts in the obvious way. It is shown in []l0[ (following the 
proof of the analogous fact for in fllf] ) that M. k is a real analytic set of dimension k(n + 1); 
on the other hand, dim X k — k(2n + 1). There are actually another k parameters, arising from the 
translations (what we called T earlier) at eachpj. Thus the correct parameter space is 2k(n + 1)- 
dimensional. In current work, D. Pollack shows that this extended version of $ is a Lagrangian 
immersion which is equi variant under the obvious action of the conformal group on either side. 

Our purpose here is to use the refined asymptotics theorem to show that the image of M k in 
X k lies in the zero set of a collection of real analytic equations obtained by computing the global 
Pohozaev integral for each solution u with respect to any conformal Killing vector field X on S n . 
Although there are not enough equations here in general to fully determine the moduli space M k , 
the relations we can obtain here are still interesting in their own right. Before embarking on this 
computation, though, we first define these integral invariants. 



6.1 Pohozaev invariants 

We now turn to a discussion of the existence and specific form of a family of homological integral 
invariants of solutions of equation (]]]). These 'Pohozaev invariants' were discovered in their 
simplest form by Pohozaev, and applied by him to prove the nonexistence of nonnegative smooth 
solutions of (jl]) satisfying Dirichlet conditions on smooth, star-shaped domains. The identity he 
discovered was later put into a natural and general Riemannian setting in |]l7| . These extended 
invariants associate to any metric g of constant scalar curvature on a manifold f2, any closed 
hypersurface S C O and any conformal Killing field X, a real number V(X, E) which depends 
only on the the homology class of E in iJ™ _1 (r2). The invariant associated to X, g and E is 

n — 2 f 

T(X, g, E) = / T(X, v) da, (22) 

2(n-l) 7 S 

where T(-, •) is the trace-free Ricci tensor for the metric g, v is the unit normal to S and da is 
the volume form induced by g on E. 



18 



In this subsection we shall compute these invariants in our particular case of interest, where 
g = u~™^5, u is a solution of (0) in the ball B(0, 1) with an isolated singularity at and £ = Ei„i 
is the sphere |z| = 7] < 1. To do this, we need two ingredients. First, the trace-free Ricci tensor T 
transforms fairly simply under conformal changes: if g = 4>~ 2 go are any two conformally related 
metrics, with trace- free Ricci tensors T and To, then 



T = T + (n - 2)<f ) - 1 ^Dd0 - ^3o) 



Next, the conformal Killing vector fields on the sphere constitute a vector space of dimension 
(n + l)(n + 2)/2. A spanning set is given by the following four basic types of conformal Killing 
fields: 

• X^ = bid Xi , the generators of parabolic motions fixing infinity, 

• T> — Y] Xjd Xi , the generator of dilation, 

• Ry^ —^2((b- x)ci — (c • x)bi) d Xi , the generators of rotations, and 

• y( b ) = ((6 • x)xi — ^\x\ 2 bi) d Xi , the generators of parabolic motions fixing zero. 

We shall compute relative to the Euclidean metric S. The transformation rule for the trace-free 
Ricci tensor yields 

T= -^—u- 2 du-du-2u- 1 Ddu- ( — ^— m~ 2 |Vu| 2 + ^-^u^ J S. (23) 
Next, the hypersurfaces T, v have unit normal and volume form 

-2 X -2 2(7i-l) , 

v = u"- 2 — = u"- 2 \x\V, and da — u "- 2 \x\ n d6. 
\x\ 

Thus we need to evaluate 

^-\\ T(X,v)da= f \ x \^l^-(X-Wu)(V-Wu)-^u(Ddu)(X,V) 
2(n - 1) J\ x \=r, J\x\=v [n-1 n-l 



1 | Vw | 2 + ^ 2)2 u ^ ) (X ■ V) \ dO. (24) 



n-l 4(n - 1) 

Denote by S the symmetric 2-tensor appearing on the right here (and including the factor l^l™ -2 ). 

The main considerations in the computations below involve homogeneity. Let us say that a 
term is of order j if it is a sum of products of terms homogeneous of order j and functions periodic 
in log Thus, from the refined asymptotics of u we see that 

S = S(_2) + <?(_x) + S', 

where Stj\, j = —2,-1, is a symmetric bilinear form, the coefficients of which (in a standard 
Euclidean basis) are of order j. All coefficients in the remainder term S' are bounded. Each of 
the conformal fields X listed above is homogeneous, and so the corresponding term Su)(X,T>) is 
of order j + i + 1 if the coefficients of X are homogeneous of order £. Since the integral itself is 
independent of the radius 77, the invariant in each case, coincides with the integral of the term of 
order zero. 

The easiest case is when X — Y^ b \ since it has coefficients homogeneous of order 2, and so 
S{X,V) = 0{\x\). We conclude that 

V(Y^ b \g)=0. (25) 
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The other cases require more work. First we make some general calculations. Using the refined 
asymptotics of u always, we first see that 



n — 2 

V ■ Vu = Ix]^ \ -(v' e H — v e ) + (a- x) 



n-2 n(n-2) 2±5 



(26) 



Next, by definition of the Hessian of a function, 

Ddu(X,V) = X{Vu) - (V x X>)u. 
Since V^D = X, this reduces to X(T> ■ Vu — u). Thus we will need that 



, 1-r. 



^ ^ - i / / n — 2. n(n — 2) , , £±f 

X> • Vu - « = |a;|— <^ -« + — — v e ) - v - ; (a • x)^- 2 + . . A . (27) 



Finally 

2n , . , , n — 2 



= M" n (vr 2 + ■ x)(-< + — Ve ) + . . .J . (28) 

First consider X = P. Since (I? ■ Vu) 2 = |x| 2 |Vu| 2 , the first and third terms combine to give 

(V ■ Vu) 2 = \x\ 2 - n ( « + ^u e ) 2 + (a-x)( «) 2 - (n ~ 2)2 + 



1 2 V -'(t/ g + — « g )j + ■■■) ■ (29) 

Combining this with the appropriate multiples of V applied to ( p7| ) and |a;| 2 multiplied by (p8f), 
we find that 

S(T>,V) = H(e)(l + n(a-x) + .. .). 
Integrating on S n , and taking the limit as i] — > 0, we get 

V(D,g)=u n _ 1 H(e). (30) 

Next, consider the case X = i?( fc ' c ) (or simply i?, for short). The third and fourth terms of 
the integrand involve the inner product of R with T>, and hence vanish. Also, R annihilates any 
function of Thus, the components of order zero in each of the remaining terms (R- Vu)(D- Vu) 
and R(T> ■ Vu — u) are of the form F(|a;|)(e • x) for some vectors e. These integrate to zero, and so 

T(R (b ' c) ,g) = 0. (31) 

For the final case, when X = X^ b \ we decompose X into the orthogonal sum of two vectors 
XO) +X ( - 2 \ where 



x 



2 



is the radial component. Then 



S{X^\V) = h -^S{V,V). 



Using the expansion above, this becomes 



b ■ x (a - x)(b ■ x) TT . . 
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The first term here integrates to zero, while for the second we use the identity 

/ (a ■ x)(b ■ x) d6 = —u> n _i(a ■ b), 

J\x\ = l n 

to conclude that its integral equals 

u) n -i(a ■ b)H(e). 

For the final term, S{X^ 2 \V), we continue with the same methods, noting that since X^ is 
orthogonal to the radial direction, the third and fourth terms once again vanish. A straightforward 
calculation, using the same sorts of parity considerations, shows that the integral of the remaining 
two terms is exactly the same as for the first component. Putting these together, we get, at last, 
that 

P(X<f'\g) = -^-u n - 1 H(e)(a-b). (32) 
n — 2 

Notice that one consequence of these calculations is that for any solution u of ([!]) with an 
isolated singularity at the origin, we can recover the values of the parameters e and a in its refined 
expansion from the Pohozaev invariants. 



6.2 The balancing formulae 

We now put the information from the last subsection together as follows. Let u correspond to 
any element of Mk- We regard u as a function on M™ rather than on the sphere S n for simplicity. 
Assume that u is singular at the points {p\, . . . ,Pk}, and has Fowler and translation parameters 
at pj given by ej , Tj and aj , respectively. Let X be any one of the conformal Killing vector fields 
listed earlier. We compute the invariants for the homologically trivial hypersurfaces 

£, = u£ =1 fl(p i ,»7)UB(0,l/t J ) 

by letting rj tend to zero. To do this, first note that u = 0(|x| 2_Tl ) as |a:| tends to infinity (since 
oo corresponds to a regular point of u on the sphere), and thus the integral around the sphere 
of radius 1/r] tends to zero with rj. As for the integrals around the other components of S, ; , we 
decompose X at each pj as a sum of these four types of basis vector fields translated by pj . Thus 

. xw = J2bid Xi 

• D = (xi - Pj,i)d Xi + PjA t 

. #cm =j2((b- Pj)ci - (c • Pj )bi)d Xi +E(0-0- pj))a -(c-(x- Pj)h)d Xi 

• y {b) = E ((b-Pi)Pa - \\P^) +E ((&■ (x- Pj )) Pji - ( Pj ■ (x-p^b^+Ub-PjK^ - 
Pji)d Xi + 0(\x -p 3 \ 2 ). 

Here pji is the coordinate of the point pj, and the last term in the last item is some vector 
field with coefficients vanishing quadratically at pj. 

Now apply the computations of the last subsection to evaluate the limit as r\ tends to zero of 
the integral around the sphere \x — pj\ = r\. For X = X^ we obtain E H(ej)(b ■ aj) = and 
since this holds for any vector b € M™, we conclude that 

k 

Y f H{e j )a J = 0. (33) 

Next, with X = V we get 

k 

Y j H(e ] ){{P J -a J ) + \) = Q- (34) 

3=1 
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Next, with X = R (b ^ we get 



k 

]T H( Sj ) ((b ■ Pj )(c ■ aj ) - (c- Pj )(b -a 3 )) = (35) 

for any vectors &, c € W 1 . Finally, with X = we get 

Y, H ^){{b-P ] )(a ] -p J + \)-\\ P A 2 {b-a ] ) \ =0 (36) 

for any 6 6 K™. These are the global balancing formula? which hold for the parameters associated 
to any solution u 6 Mk- 

This set of analytic equations does not determine all the parameters of the solution u, and in 
particular, gives absolutely no information about the 'translation parameter' Tj appearing in the 
asymptotics formula at each pj. Still, particularly when k is small relative to n, these equations 
do shed more light on the global nature of the moduli space. For example, when k = 3 this set 
of equations can be solved and leads to an explicit formula for the parameters aj in terms of the 
£j and pj. This is not possible in general, but nevertheless we still have the following result: 

Proposition 6 There exists a constant C > 0, depending only on A = {pi, . . . ,Pk}, such that 
for any u € .Ma, the corresponding translation parameters and necksizes satisfy \H(sj)aj\ < C , 
j = 1, . . . ,k. 

We have already seen that V{X^ h \g, dB(pj,n)) determines w„_i(aj -b)H(ej). On the other hand, 
using the universal upper bound of Theorem |2| and elliptic estimates, we also get an a priori bound 
for this invariant. Since this is true for all b £ K™, the result follows. 



7 Nondegeneracy of the moduli spaces near their ends 

In this final section we give another application of the asymptotics theorem and address the issue 
of the nondegeneracy of the unmarked moduli space Mk- Unlike in the previous sections, we 
use only the simpler asymptotics result, not the more refined one. We have included this here 
because the arguments are soft, and not too different in spirit from some of the ones used above. 
This result is an adaptation of one in Q , and our desire here is to show its validity beyond the 
more limited setting of that paper. 

It is unknown whether degenerate solutions ever exist. If they do they are quite unstable: it is 
shown in pd[ | that under an arbitrarily small generic change in the conformal class [go] , the moduli 
space becomes smooth. On the other hand, it is also nontrivial to show that a given solution 
is nondegenerate. This should be easier with explicitly constructed solutions, but unfortunately 
this was still impossible to do with the first-known solutions from Jl7|| . One construction of 
nondegenerate solutions, in a somewhat limited setting, was given in ]12|, and another much 
more general one was given in |J. In this last paper it was shown that given any singular set 
A, there is a nondegenerate solution singular at the points of A; these solutions have very small 
necksizes (i.e. Fowler parameters). (It was also shown that for generic configurations A, these 
solutions are also nondegenerate in the marked moduli space .Ma-) 

We show here that an argument from ^| may be adapted to prove something slightly weaker, 
although probably optimal. Before we state it, we discuss briefly the compactification theory of 
these moduli spaces. It was shown in jn| that if ug is any sequence of elements in A4k such that 
the singular points p^ are bounded away from one another and all Fowler parameters e!- are 
bounded away from zero, then some subsequence of ug converges to an element Uoo € Mk- This 
result limits the ways in which noncompactness in Mk can occur. This was discussed further 
in and pOL where it was shown that if ug is a sequence in Mk with one or both of these 
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restrictions not satisfied, then it is possible find conformal transformations Fg such that some 
subsequence of F^ug converges to an element u 1 ^ in some My with k! < k, or else converges 
to zero uniformly on compact sets. (To make things consistent here, we let Mo denote the set 
of pullbacks of the standard (smooth) metric on S n by conformal transformations.) This result 
states, then, that Mk may be compactified by adding to its ends certain subsets of moduli spaces 
Mk> of solutions with fewer or no singular points. 
Finally then we can state our result. 

Proposition 7 Let ug G Mk be any sequence of elements such that the singular points stay 
bounded away from one another. Suppose that this sequence converges to G Mk' with k' < k. 
Then either Uoo is a degenerate solution in Mk' or else ug is nondegenerate for sufficiently large 
I. 

For simplicity in the notations, we will assume that the singularities {p^\ ■ ■ ■ >pi } do not depend 
on I. 

The proof is by contradiction. First let us apply the refined asymptotics theorem as follows. 
Choose small balls B(pj, p) which are disjoint from one another and such that we may write each 
U£ as a sum of two functions 

§,a,T . 

ug = u£ ' + we. 

Here e = {e[ e) , ...,e^},T = {T^\ T^} and a = {a[ e \ . . . , a^} are the Fowler and trans- 
lation parameters at each pj for ug and u e l ' aT is a function agreeing with ug outside the balls 
B(pj,p) and equalling the model deformed Fowler solution (relative to the background spherical 
metric) u e ^ a ..Tj in B(pj, \p); finally, wg G C~ a for some fixed 7 > (4 — n)/2 close to (4 — n)/2. 
Assume that at least some of the ey^ tend to zero; otherwise the theorem is trivial. Relabel the 

points so that, after passing to a subsequence, Ej tends to zero for k' + 1 < j < k while ey 
converges to some nonzero values Ej for j < k' . Finally, assume that there exists, for each £, a 
function (fig G C?' Q such that he(f>e = 0, where is the Jacobi operator at ug. 

Although it is not literally true, we shall assume that the singular points pj do not vary with 
£. Since we are assuming that they stay a bounded distance away from one another, we could 
transform to this case where A is fixed by a convergent set of diffeomorphisms of the sphere, but 
the only effect this would have would be to complicate notation. 

Normalize <pg, multiplying it by a suitable constant, so that sup d(y)~ 7 \ (j)e(y)\ = 1, where d(y) 
is the distance of the point y from the singular set A in the spherical metric. Choose a point 
ye G S n \ A realizing this supremum, i.e. such that d{ye)~ 1 \4>e{ye)\ = 1< 

If some subsequence of the ye converges to a point yo € S n \ A, then we may extract a 
subsequence of the ve converging to an element v' G S n \ A', where A' = {p\, . . . ,pk'} and also 
so that (fie converges to a nontrivial function <fi on S n \ A'. Clearly, \<fi\ < d{y)^ , and also U<fi = 
weakly on all of S n , where L' is the Jacobi operator at v' . Since 7 > (4 — N)/2 and v' is smooth 
at the points Pj, j > k' ', it follows from a standard removable singularities theorem that <fi is 
actually smooth across these points. Thus <f> G C~' a (S n \ A'), and so we have shown that the 
limiting solution v' G Mk', is degenerate. 

If, on the other hand, some subsequence of the ye converges to one of the points pj , then it is 
more convenient to transform the problem, using its conformal equivariance, to one on a cylinder 
before proceeding further. First, choose a function A on S n \ {pj,q} (for any second point q) 
such that A~^g is the product metric gc = ^^(dt 2 + dO 2 ) on the cylinder C = lx S n ~ x , 
with t = 00 corresponding to Pj and t = —00 corresponding to q. On C, the function A is simply 

a multiple of (coshi) - 5~, and on S n is of the order dist(y, {pj, q})~ *~ ■ The solutions vg on S n 
correspond to solutions Ave on C . Since the metrics gc and go both have scalar curvature n(n— 1), 
a straightforward calculation (cf. Jl^|) shows that the linearized scalar curvature operators hg on 
S n at vg and hc,e on C at Ave satisfy the same conformal equivariance property as the conformal 
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Laplacians for these two metrics, namely 

JV+2 

Lc,i(A(j>) = i^L^, (37) 

for any function <p. (This depends strongly on the fact that both metrics have the same scalar 
curvature.) 

Using this transformation, we now have a sequence Ave of solutions of the nonlinear equation 
as well as a sequence Acj>e of solutions of the Jacobi operator at Ave on C. For simplicity we 
relabel these functions and the operator by ve, <fii an d again. Let y = (t, 8) denote the variable 
on the cylinder and define 7 = 7 + (n — 2)/2, so that 7 > 1, close to 1. Then 

sup%)-^| = l, (38) 

where d(y) is once again a smoothed distance function to the singular points in A \ {pj}, trans- 
planted to C, in some large compact set and equal to sechi outside this neighbourhood. By (pq), 
4>£ decays at both ends of the cylinder. 

As before, let the supremum in ( |3"s| ) be attained at the point ye — (te,9e)- By assumption, 
te — * 00. Translating back by te and renormalizing the solution, we find yet another sequence 
of solutions, which we again call we, attaining their maximum at t = 0, and which solve the 
translated equation, which we again write as \*i<\>e — 0. Here he is the Jacobi operator at 
ve(t + te,0). As usual, some subsequence of the <pe converge to a nontrivial solution <f> of the 
limiting equation L(j) = 0, and <f> is bounded by e~ 7 * for all t. 

There are two cases to consider. In the first, pj is one of the singular points for which s), 
tends to zero. Here there are two subcases, depending on whether ve(te,0e) tends to zero or not. 
If it does tend to zero, then <j) satisfies the equation 

n /„n . \ n(n — 2) , n ( ' . (n — 2) 2 \ , 

Decomposing <fi into its Ag eigencomponcnts, we see that any eigencomponent cj>j is a sum of 
exponentials, <j>j — a^e^ jt + aje^^ ljt . Since <j> decays as t — > +00, = 0. But then it is clear 
that no function of the form e _Mj t can be bounded for all t by e -7 ' unless 7 = fij, which is not 
the case, so we arrive at a contradiction. In the other subcase, ve(te,9e) does not tend to zero. 
Translating by a fixed finite amount, we may assume that ve tends to the function (coshi)^"^/ 2 , 
and hence, after pulling out the superfluous constants, that the limiting function (f> satisfies 

( „ 2 a (n-2) 2 n 2 ~A l2 \, 

f d 2 t + A e - y ' + ^^sech 2 n = 0. 

Again separate cj> into its eigencomponents 4>j- Then 

/ (n — 2) 2 \ n 2 — 4 
d 2 th ~ {^-^^ + ^ + -^sech 2 i ^ = 0. 

For j = the indicial roots of this equation at both ±00 are ±(n — 2)/2, for j = 1 they are ±n/2, 
and for j > 1 the indicial roots are all > [n + 2)/2. 

The components 4>j with j > 1 are easy to eliminate. In fact, these 4>j must decay faster than 
e ±(«+2)|t|/2 _|- 0Oj so we ma y multiply the equation satisfied by (f>j and integrate by parts to 
obtain 

J {dtWj f + (\j + 4 J 7— sech 2 t (j) 2 dt = 0. 

Since Xj > 2n, the integrand is nonnegative, hence <j)j = 0. 

For the remaining cases, when j = 0,1, the indicial roots at ±00 are less than (n + 2)/2 
in absolute value. On the other hand, to check unmarked nondegeneracy it suffices to use any 
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7e 3 ,n+i > 7 > 1. But because ej tends to zero, this upper limit tends to (n + 2)/2. Thus if 
we choose 7 in the range (n/2, (n + 2)/2) then both 4>o an d <t>i decay less quickly than e~ 7 * as 
t — ► +00, which implies that these (f)j too must vanish. This is a contradiction. 

The final case to consider is when Sj does not converge to zero, so that vg is converging 
to some Fowler solution v s . In this case 4>g converges to a solution of C e 4> = which satisfies 
\4>\ < Ce~ 7t . But we have shown that all solutions of this equation are sums of terms for each 
eigencomponent of Ag which satisfy bounds |e =F7e j * < C. Clearly this is incompatible with 
the previous bound, so we arrive here too at a contradiction. 

We have shown, finally, that unless it is converging to a degenerate solution in some Mk'y v t 
is nondegenerate in M.k- 

It may well seem disappointing that we can not exclude degeneracy near any end, but as 
remarked earlier, this form of the result is probably optimal. However, in certain cases we can 
deduce nondegeneracy without restriction; this is due to the fact that only the first of the several 
cases treated in the proof did not necessarily lead to a contradiction. However, for example, if 
the solutions ve are converging to zero uniformly on compact sets of S n \ A, then they must be 
nondegenerate for sufficiently large £. This is because the equation satisfied by the limiting Jacobi 
field (f> is (Ag™ — (n — 2) 2 /4)0 = which has no nontrivial solution. This is the case studied in 
. Another case where we can deduce nondegeneracy for sufficiently large I is when all but two 
of the singular points pj disappear in the limit, i.e. when k' — 2. This is because any Fowler 
solution is nondegenerate. On the other hand, it is not possible, using this argument, to deduce 
nondegeracy for the solutions constructed in Jl7[ | , for there vi converges to the constant function 
1 on S n , and the sphere is degenerate. 
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